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Abstract. We have analysed the angular power spectra of the Parkes radio continuum and polarisation survey of 
the Southern galactic plane at 2.4 GHz. We have found that in the multipole range / = 40-250 the angular power 
spectrum of the polarised intensity is well described by a power-law spectrum with fitted spectral index oil = 
2.37 ± 0.21. In the same multipole range the angular power spectra of the E and B components of the polarised 
signal are significantly flatter, with fitted spectral indices respectively of ole = 1-57 ± 0.12 and olb = 1-45 ± 0.12. 
Temperature fluctuations in the E and B components are mostly determined by variations in polarisation angle. 
We have combined these results with other data from available radio surveys in order to produce a full-sky toy 
model of Galactic synchrotron intensity and linear polarisation at high frequencies (y <; 10 GHz). This can 
be used to study the feasibility of measuring the Cosmic Microwave Background polarisation with forthcoming 
experiments and satellite missions. 
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1. Introduction 

Current cosmological models predict the level of polarisa- 
tion of the Cosmic Microwave Background (CMB) to be 
at 5-10% of the temperature anisotropies at the 20 arcmin 



dial density fluctuations and on the thermal history of 



angular scale, that is at fiK level (e.g. Bond & Efstathiou 



1987). These temperature variations have to be detected 



over a 3 K continuum with instrumentation whose typ- 
ical noise temperature is of the order of tens of Kelvin. 
The detection of the CMB polarisation represents there- 
fore a considerable challenge. To date, only upper limits on 
the CMB polarisation exist. The most recent results from 
ground based observations place an upper limit of 10 /iK 
on any polarised signal from the sky at the 7° angular 



the universe (e.g. Kosowsky 1996). The variation of po- 
larised CMB emission is predicted to peak at sub-degree 
angular scale. The NASA satellite mission MAP will have 
enough sensitivity to detec t the CMB polarisation, if this 
is present at the uK level ( Kogut fe Hinshaw 2000] ). The 
future ESA mission Planck will have higher sensitivity and 
may provide direct imaging of the CMB polarised signal 
flMandolcsi fc Puget 1998)) . The future SPORT experi- 
ment that is planned to be placed on the International 
Space Station also has the sensitivity to detect /iK-level 
polarised sig nal from the CMB if this is present at angular 
scales Z 7° ( |Fabbri et al. 1999| ). 



scale in the frequency band 26-36 GHz ( Keating et al 



2001 ) and an up per limit of 15 /iK at the l°-0°.5 angular 
scale at 90 GHz ( |Hcdman et al. 200l[ ). 

Current cosmological models also predict a shape for 
the angular power spectra of the polarised component of 
the CMB, so the measure of these spectra will provide 
additional information on the properties of the primor- 
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A fundamental question for any measurement of the 
polarisation of the CMB is whether the expected signal 
can be distinguished from the foreground polarised signal 
from our galaxy. The interstellar magnetic held is illumi- 
nated by cosmic-ray electrons which spiral around held 
lines and thereby emit synchrotron radiation. This radia- 
tion is intrinsically highly lin early polarised, 70-75% in a 
completely regular field (e.g. Rybicki fc Lightman 1979). 



Thermal emission from dust may also be highly polarised, 
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depending on the shape and alignment of the dust parti- 
cles flWright 1987[ |Prunet et al. 199S| ). 

As in the case of the unpolarised emission, knowledge 
of the angular power spectra of the foreground polarised 
component is essential when trying to recover the an- 
gular power spectra of the CMB. In this paper we use 
the Parkes 2.4 GHz polarimetric survey of the southern 
Galactic plane in order to derive the angular power spec- 
tra of the polarised component of the Galactic synchrotron 
emission. The survey has relatively small sky coverage but 
has 10 arcmin angular resolution, which allows the statis- 
tical properties of polarised synchrotron emission to be in- 
vestigated at the angular scales where the CMB polarised 
signal is expected to peak. 

In order to assess the feasibility of CMB polarisation 
measurements, study the effect of instrumental systemat- 
ics and devise the data analysis strategy it is necessary 
to perform simulations of the satellite missions. For these 
simulations an input sky with all the known components 
of the emission is necessary. By exploiting existing obser- 
vations of the microwave sky, we have to strive to con- 
struct a simulated sky which represents our best guess of 
some of the properties of the real sky in the unexplored 
regions and frequency channels. This will allow us to per- 
form more realistic simulations of the mission and to be 
in a better position to interpret the real data when these 
become available. 

We use the results of our analysis of the Parkes survey, 
combined with other data on synchrotron emission coming 
from total intensity large sky radio surveys, to produce a 
full-sky synthetic map of synchrotron polarised emission. 
The paper is organised as follows. In Sect. || the global 
angular power spectra of the Parkes 2.4 GHz polarimetric 
survey are derived. The simulations performed to assess 
the reliability of the spectra are presented in Sect. ||. The 
results are discussed in Sect. f|. In Sect. ||the derived spec- 
tra are used as basis to construct a toy model of the all-sky 
polarised synchrotron emission. 

2. The Parkes data at 2.4 GHz 

The Parkes 2.4 GHz survey is a sensitive, polarime tric sur- 
vey of the southern Galactic plane carried out by Duncar| 
|et al. (1997)| . The survey covers 127° of Galactic longitude 
(238° < I <5°), with a latitude coverage out to at least 
±5° (up to b =7° and b =8° over some longitudes). The 
angular resolution of the images is 10.4 arcmin. The nom- 
inal rms noise of the survey (for the total and polarised 
intensities) is 5.3 mK. The survey is publicly available and 
the data come split into 6 rectangular fields sampled at a 
pixel resolution of 4 arcmin. The linear polarisation signal 
is given in terms of the two Stokes parameters Q and U. 

The Parkes survey has been recently anal ysed by 
Baccigalupi et al. (2001)1 and |Tucci et al. (2000)| . In both 
of these previous studies, the data were split into twelve 
square patches, each of size 10°xl0°, which were then in- 
dividually analysed. Our work differs from these in that 
the power spectra of the different polarisation components 



are derived from a global analysis of the complete Parkes 
data. 

In order to perform this global statistical analysis of 
the survey we have resampled the six fields of the original 
data into HEALPix tessellations. The iso-Latitude prop- 
erty of the HEALPix tessellation of the sphere makes it 
ideal for fast computation of the angular power spectrum 
of a field defined on the sphere ( Gorski et al. 1999| ). In 



the HEALPix tessellation the angular size of the pixels is 
determined by the nside parameter. We used a HEALPix 
tessellation with nside=102A which corresponds to pixels 
of linear size 3.4 arcmin. When resampling the survey, the 
pixel value in the HEALPix tessellation was assigned by 
weighted integration of the values of the original pixels 
that overlapped with the new pixel, where the weight of 
each pixel contribution is the fraction of its area inter- 
secting with the new pixel. This approach ensures that no 
spurious pixel-to-pixel signal correlation is introduced by 
the resampling. 

The three channels of the survey (T, Q and U), resam- 
pled into HEALPix are shown in Fig. p\. Th e derived inten- 
sity of linear polaristion, L = y/ Q 2 + U 2 , is also shown. 
The maps are in Galactic coordinates. 

The most striking difference between the total inten- 
sity channel (T) and the polarisation channels (Q, U and 
L) is that the total intensity signal is dominated by the 
signal from the Galactic plane while in the polarisation 
channels the Galactic plane is hardly visible. Moreover 
many point sources which are clearly visible in the total 
intensity channel do not exhibit a polarised counterpart. 

In the polarised intensity image the most apparent fea- 
tures are the extended areas of polarisation, with sizes of 
the order of 5° or so, which do not appear to be con- 
nected with discrete sources of total power emission. Such 
regions lie around longitudes 265°, 280°, 320° to 330° and 
355° to 5°. While one of these structures can be iden- 
tified with the Vela supernova remnant at I = 265°, the 
origin and nature of the other features is not yet explained 
( puncan et al. 1997| ). 

Five regions that are bright in the T maps have been 
blanked out to prevent a significant fraction of instrumen- 
tal polarisation from being introduced. These are Sgr A 
and the bright Hn complexes G353.2+0.8, G284.14-0.3, 
G287.5-0.6 and G291.28-0.71. 



2.1. Derivation of the power spectra 

The three Stokes parameters T, Q and U fully describe 
any state of linearly polarised light. However, whereas the 
temperature is a scalar quantity (invariant under rotation 
of the plane of the sky), Q and U are not. They depend 
on the direction of observation n and on the Cartesian 
axes perpendicular to h used to define them. Recent theo- 
retical developments have shown that it is more useful to 
describe the polarisation field in te rms of two quantities 
that arc invariant under rotation ( Zaldarriaga fc Seljak 
1997 ; Kamionkowski 1997 ). These two quantities are usu- 
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Fig. 1. The Parkes survey at 2400 MHz. The six fields of the original data have been resampled into one HEALPix 
tessellation with a pixel size of 3.4 arcmin. The first three images (from the top) give the three Stokes parameters of the 
signal: T, Q, U, respectively; the lowest image shows the derived polarised intensity L = \J 'Q 2 + U 2 . The temperature 
scales for T and L have been histogram equalised. 



ally called E and B and they are obtained by a linear, 
non-local, transformation of the Q and U map (hereafter 
we will refer to this transformation as to the E-B trans- 
form). The E and B components differ in their behavior 
under parity transformation: B changes sign while E does 
not. Four power spectra are needed to characterise the 
fluctuations in a Gaussian theory: the autocorrelation of 
each of T, E and B and the cross correlation between E 
and T. 

Current cosmological models predict the CMB inten- 
sity and polarisation variations to be a Gaussian field and 
therefore give predictions for the four power spectra T, E, 
B and E xT. Galactic diffuse emission is not a Gaussian 
field and so cannot be fully characterised by the four power 
spectra, however, since we are interested in comparing the 
statistical properties of this emission with predictions for 
the CMB power spectra we have derived the four power 
spectra of T, E, B and E x T of the Parkes data. 

The power spectra have been obtained by perform- 
ing a spin-weighted harmonic analysis of the field in the 
HEALPix grid with the ANAFAST program (the soft- 
ware is part of the HEALPix package) . The results of the 
analysis for the whole set of data are summarised in the 



following subsections. The angular power spectra of the 
polarised intensity was derived by simple spherical har- 
monic decomposition, also using the ANAFAST program. 

From the analysis of large sky surveys, the angular 
power spectra of the various components of unpolarised 
Galactic foreground emission appears to be well described 
by power law spectra with a spectral index between 2 and 
3. Below, we model the angular power spectra of the dif- 
ferent components of the Parkes data with power-laws of 
the form: 

Cf = kl~ ax 

where I is the multipole order and X may stand for T, L, 
E, B, e and 6, with e and b indicating the E-B transform 
of the sinus and cosinus of the polarisation angle. 

The limited sky coverage of the survey and the FWHM 
resolution of 10.4 arcmin limit our spectral analysis to the 
multipole range / = 40-250. Sect. || describes Monte Carlo 
simulations that were made to assess the effects of reduced 
sky coverage and finite resolution. 
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Fig. 2. The global angular power spectra of T and L = 
\/ Q 2 + U 2 for the Parkes survey (continuous lines). The 
power spectra of T and L after median filtering are also 
shown (dashed lines). Point source suppression by the fil- 
tering process is responsible for the steepening of the in- 
tensity spectrum. Such steepening is not observed in the 
spectrum of polarised emission (up to I ~ 150). The hor- 
izontal line gives the nominal noise level. The multipole 
order I corresponds to a typical angular scale of = it /I. 



2.1.1. Total and polarised intensity 

The angular power spectrum of the T channel of the 
Parkes survey is shown in Fig. 0, In the multipole range 
I =40-250 a linear least-squares fit to this spectrum gives: 
a?T = 1.67 ± 0.150. In the figure the angular power spec- 
trum of the intensity of the polarised signal L is also 
shown. Linear least-squares fit to this spectrum gives 
a L = 2.37±0.21. 

From the images in Fig. |l| point sources are clearly 
identifiable in the map of total emission while they are 
not identifiable in the polarised channels. In order to ver- 
ify whether point sources or noise at pixel scale is affecting 
the spectra derived from the Parkes survey, we median 
filtered the intensity map (T) and the derived polarised 
intensity map (L) by convolving the data with a kernel 
of 9x9 pixels. We have used this technique to derive the 
angular power spectrum of diffuse emission from the large 
sky radio survey of lonas et al. (1998) . We showed that 
median filtering can be used to effectively remove point 
sources and derive the angular power spectrum of the dif- 
fuse component (Giardino et al. 2001). 



The angular power spectrum of a population of point 
sources randomly distributed in the sky is Poissonian, that 
is with constant power at all I. Therefore, if point sources 
contribute significantly to the temperature fluctuation in 
the field of sky, the power spectrum of that field is ex- 
pected to steepen when the point sources are removed and 
only the diffuse components remain. The same happens if 
the signal is dominated by instrumental random noise. On 
the other hand if the signal in the field is dominated by 
the diffuse component, median filtering will not have any 
effects on the spectrum of the emission up to the mul- 
tipole order where the high-frequency fluctuations of the 
emission are averaged out. We estimated by Monte Carlo 
simulation that, for a beam with a FWHM of 10.4 arcmin 
and a map with pixel size of 3.4 arcmin, median filtering 
with a 9x9 pixel kernel allows the input angular power 
spectrum of a diffuse component (with a spectral index of 
a = 2) to be derived up to I — 150 with an error lower 
than 3%. 

The results of the analysis of the median filtered T and 
L maps are also shown in Fig. ^|. 

From the figure it is apparent that the spectrum of the 
total emission steepens significantly upon the application 
of the filtering process. The fitted spectral index values, in 
the / range 40-150, before and after median filtering are 
respectively a/ = 1.68 ± 0.03 and a} nf = 2.05 ± 0.03 (only 
the formal fit errors are given in this caseQ). This indicates 
that a significant fraction of the power spectrum at inter- 
mediate and high I is due to the presence of discrete sig- 
nals: either a population of point sources or instrumental 
noise at the pixel scale. In fact instrumental noise at pixel 
scale is unlikely as the nominal noise level is well below 
the spectrum derived for the total intensity. We therefore 
interpret the steepening observed as due to the removal of 
the point source contribution. 

Fig. U shows no comparable change in the slope (for 
I < 150) of the angular power spectrum of the derived 
L map after median filtering. The fitted spectral index 
values before and after median filter are respectively at = 
2.34 ± 0.03 and = 2.45 ± 0.03. 

This analysis implies that, up to multipole order I = 
150, the fluctuations in the polarised emission are of dif- 
fuse nature and they are not due to the presence of discrete 
signals such us point sources or instrumental noise at pixel 
scale. 

Even though we can exclude contamination of the po- 
larisation spectrum by point sources and instrumental 
noise at pixel scale only up to multipole I =150, in the fol- 
lowing sections we derive the spectral indices for the other 
components of the Parkes data over the entire multipole 
range I = 40-250 that can be directly probed given the 
survey's sky coverage and angular resolution. The reader 
should note from Fig. || that the slope of the power spec- 
trum of the unfiltered polarised emission does not obvi- 



1 Errors quoted in this paper are estimated via Monte Carlo 
simulations unless otherwise stated 



2 Since in this case we are interested in the spectral indices' 
relative values before and after the filtering, the formal fit er- 
rors are representative of the uncertainty in this comparison 



G. Giardino et al.: Towards a model 



5 



ously change between an I of 150 and 250, and therefore 
the fits given up to 250 should also be consistent with 
those made only up to 150. Indeed, we have explicitly es- 
tablished this. 

2.1.2. The E and B components 



6[arcmin] 
1000 100 10 




10 100 1000 

1 [multipolej 



Fig. 3. The power spectra of the E and B components 
for the Parkes data. The linear least squares fits to the 
spectra in the multipole range 40-250 are also shown. 



The E x T spectra derived from the ten simulations 
showed a standard deviation around of ct; = 1.4 • 
10~ 4 /Z(7 + 1) K 2 . The grey shaded region in Fig. [| shows 
the area within ±3crj. The comparison of the data points 
with this area is consistent with the absence of correlation 
between E and T. If a level of correlation between E and 
T exists, this is below the level that can be probed with 
the Parkes data, given the survey's sky coverage and beam 
FWHM. 

The power spectrum derived from E and the median 
filtered T is also shown in the figure. The level of correla- 
tion is hardly affected by the filtering process (up to the 
multipole order Z = 150 where the two spectra are directly 
comparable), implying that the lack of correlation between 
the T and the E component is not due to the presence of 
unpolarised point sources (which only affect the T map). 



0[arcmin] 




E xT 

E xT (median filtered T) 



The power spectra of the E and B components de- 
rived from the Parkes survey are shown in Fig.[|. In the 
multipole range I = 40-250 linear least-squares fits to the 
spectra give = 1.57 ± 0.12 and olb — 1-45 ± 0.12 re- 
spectively for E and B. As can be seen from the figure the 
power spectra of the E and B components are very similar 
in slope and amplitude. This is expected if the observed 
emission field does not show clear patterns or symmetries 
(|Zaldarriaga 200 1|). 



-0.001 — 

50 100 150 200 250 300 

1 [multipole] 

Fig. 4. The cross-correlation E x T for the Parkes data. 
The grey area gives ±3tr from zero of the recovered E xT 
spectra for 10 simulations with input null ExT (see text). 



2.1.3. The ExT component 

The cross-spectrum of the E and T components is shown 
in Fig. | 

In order to assess whether the computed correlation 
signal corresponds to null correlation between E and T 
we performed a series of simulations. We generated ten 
random Gaussian polarised fields using the spectral in- 
dices and normalisation derived from the Parkes data to 
define the power spectra of T, E and B and setting C / ExT 
to zero. The same sky coverage and beam FWHM of the 
Parkes survey were reproduced (see Sect. ||). 



2.1.4. The polarisation angle 

In order to quantify the statistical properties of the polari- 
sation angle variations we constructed the maps of cos(2#) 
and sin(2(9), where 9 is the polarisation angle. These maps 
are simply given by Q / 1 L and U / L respectively. The sinus 
and cosinus maps were then E-B transformed to obtain the 
rotationally invariant quantities, e and b. The power spec- 
tra of e and b were also derived by spin- weighted harmonic 
analysis. The result of the analysis is shown in Fig. |[ 

In the multipole range I = 40-250, the power spectra of 
e and b are well described by power law spectra with fitted 
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Fig. 5. The power spectrum of the E-B transform of 
cos(2#) and sin(2(9), where 9 is the polarisation angle. The 
transformed map are here indicated with e and b. The lin- 
ear least squares fits to the spectra in the multipole range 
/ = 40-250 are also shown. 



I < 40. Above multipole order I > 250 the derived power 
spectra are affected by the beam cut-off. 

In order to quantify the accuracy of the recovery we 
have computed the best-fit power law of the spectra ob- 
tained from the f realisations, in the two different cases of 
sky coverage. The spectral indices of these best-fit power 
laws can be compared with the true spectral indices of the 
power spectra used as inputs for the 10 realisations. The 
standard deviation of the derived spectral indices with re- 
spect to the true (input) values provides an estimate of 
the standard error on the spectral index. The results are 
summarized in Table [l]. In the case of full sky coverage 
and considering the multipole range I = 3-250, one can 
recover the input spectral indices of T, E and B with an 
error of approximately 2% (as can be seen by taking the 
ratio of the SD over the mean in three columns of the 
table) . The limited sky coverage of the Parkes survey im- 
plies that one can only fit a power law to the recovered 
spectra in the multipole range I — 40-250. By considering 
the spectrum in this limited multipole range the error on 
the recovered spectral indices increases significantly even 
when using spectra derived from full-sky maps (around 
5%, as taken from Table [l]). The error is around 10% 
when recovering the indices of the T, E and B spectra 
from maps with the limited sky coverage of the Parkes 
survey. 



spectral indices of a e = 1.74 ± 0.14 and ctb = 1.69 ± 0.13 
respectively. These values are very close to the spectral 
indices of the power spectra of the E and B components. 

3. Error assessment 

In order to verify the reliability of the power spectra de- 
rived from the Parkes survey we performed a series of 
Monte Carlo simulations. We generated 10 full-sky realisa- 
tions of a random Gaussian polarised field defined by four 
power spectra with shape and normalisation similar to the 
ones of the Parkes data. For T, E and B we used power- 
law spectra with spectral indices of ar = 1-7, Q£ = 1.5, 
olb = 1-5 respectively, for E x T the spectrum was set to 
zero. The 10 full-sky random realisations of these spectra 
were convolved with a Gaussian beam with a FWHM of 
10.4 arcmin and sampled using a HEALPix tessellation 
with nside= 1024. The power spectra for each realisation 
were then derived by performing spin-weighted harmonic 
analyses for both full-sky coverage and the limited Parkes 
sky coverage. 

The averaged power spectra from the 10 simulations 
are shown in Fig. |^. The top panels shows the derived 
average power spectra of T and the bottom panels the de- 
rived average power spectra of E and B. The two different 
cases of sky coverage are shown. The effect of the incom- 
plete sky coverage is visible in the spectra in the right 
panels where the average of the derived spectra deviates 
significantly from the input spectra (the true value) for 



4. Comparison with previous results and 
discussion 



The spectral index of ar = 1.67 ± 0.15 that we derived 
from the Parkes data for the total intensity (T) is signif- 
icantly lower than the spectral index values derived from 
large sky radio survey (Bouchct & Gispcrt 1999; Giardinc 
|ct al. 20"01 ). However this value refers to the raw data 
which beside the diffuse emission contains a significant 
population of point sources. After point source removal 
the index steepens to olt — 2.05 ± 0.15. This value can be 
compared with the value of olt = 1-9 for the spectral in- 
dex of the same region of the sky of the Rhodes/HartRAO 
survey after point source removal. 

The spectral index of the Rhodes/HartRAO survey 
steepens with Galactic latitude from a = 2.4 when all 
the data are considered to a = 2.9 when only the data 
at galactic latitude |6| > 20° are taken into account 
(Giardino et al. 2001). Therefore it is not surprising that 
when a region within Galactic latitude \b\ = 8° is consid- 
ered the derived spectral index is flatter. The reason for 
this is probably a combination of the contributions from 
supernovae remnants and from the high concentration of 
diffuse H n emission regions in the Galactic plane. 

As earlier noted, the Parkes surve y wa s previously 



analyse d by |Baccigalupi et al. (2001)| and |Tucci et al 



(2000)| . Baccigalupi et al. restricted their analysis to the 



total intensity and the polarised intensity data (L) while 
Tucci et al. (2000) present results for the power spectra of 
T, E and B, but not L. 
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Fig. 6. The average spectra of T (top) and of E and B (bottom) obtained from the spin-weighted harmonic analysis 
of 10 different realisations of a polarised Gaussian radiation field (continuous lines). The dashed lines indicates the 
input power spectra for T, E and B. These are power law with spectral indices olt = 1.7, ole = 1.5, as = 1.5 (very 
similar to the Parkes data). The input spectrum for E x T was set to zero. 



Baccigalupi et al. (2001) obtained an average value of 
at, = 1.7 ± 0.2 for the spectral index of the angular power 
spectrum of the polarised intensity of the twelve patches 
in the Z-range 100-800. This value is significantly lower 
than the value of = 2.37 ±0.21 that we derived for the 
multipole range I = 40-250. However, there are two issues 
to consider: 



1. the angular power spectrum of the unfiltered L (Fig. 
shows a flattening of the slope around I — 300. 

2. in order to analyse the Parkes data beyond I ~ 250 
where the beam suppression becomes appreciable, as- 



sumptions about the noise amplitude and distribution 
are necessary. 



If we perform a linear least squares fit to Cf 1 in the 
multipole range I =100-800, assuming purely Gaussian 
random noise and the nominal noise level given in Duncan 
et al. (1997), we obtain az, 1 i=ioo-80o = 1-75, in close agree- 
ment with the value derived by Baccigalupi et al.. On the 
other hand, if we take into account the presence of the 
slope break and fit the curve with two power laws, we 
obtain ol,z=40-250 = 2.35 and ol,z=250-800 = 1.64, indi- 



cating a flattening of the spectrum at I > 250. 
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OLT 




OLE 




OLE 




I range 


Sky coverage 


mean 


SD 


mean 


SD 


mean 


SD 


3-250 


Full 


1.70 


0.02 


1.51 


0.03 


1.51 


0.03 


40-250 


Full 


1.76 


0.08 


1.56 


0.08 


1.55 


0.08 


40-250 


Parkes 


1.78 


0.15 


1.57 


0.12 


1.55 


0.12 



Table 1. Mean and standard deviation of the recovered spectral indices from 10 different sky realisations of a random 
Gaussian field with input spectra with spectral indices ar = 1.7, Q£ = 1.5, olb — 1.5 for T, E and B. The spectrum 
of E x T was set to zero. Two different cases of sky coverage have been considered: full sky coverage and the Parkes 
survey sky coverage. 



Clearly, Baccigalupi et al. derived a flatter spectral in- 
dex by fitting the spectrum of L in the Z-range 100-800 
without accounting for the slope break. As to the rea- 
son why their analysis did not see evidence for a change 
in slope of the power spectrum of L, it is likely that by 
analysing small patches, the derived individual spectra are 
noisier, thus rendering it more difficult to detect such a 
change. 

Unfortunately, assessing the significance of this flatten- 
ing is difficult precisely because the analysis for I > 250 
is sensitive to assumptions about the noise amplitude and 
its distribution. 



Tucci et al. (2000) derived the average angular power 
spectra of the E and B components for the twelve patches 
and concluded that in the multipole range 36 < I 5s 1000 
they are both well approximated by power laws with spec- 
tral index a ~ 1.4. This value is very similar to the values 
of ole = 1.57 ±0.12 and a B = 1.45 ±0.12 that we derived 
for the E and B components respectively. 

The angular power spectrum of e and 6, the E-B trans- 
form of the sinus and cosinus of the polarisation angle, 
were first derived in the present work. In the multipole 
range I = 40-250 the power spectra of the e and b com- 
ponents of the Parkes data are well described by power 
laws with fitted spectral indices of a e = 1-74 ± 0.14 and 
Qfb = 1.69±0.13. These values are very close to the spectral 
indices derived for the E and B components. Fluctuations 
in e and b are only determined by variations in polarisation 
angle. Therefore this analysis shows that fluctuations in 
the E and B components are mostly determined by varia- 
tions in the polarisation angle rather than polarisation in- 
tensity, which, over this multipole range, is characterised 
by a significantly steeper power spectrum. 

Finally, we did not find evidence of correlation between 
E and T in the Parkes data. This lack of correlation can 
be interpreted as a direct consequence of the fact that fluc- 
tuations in the E components are mostly determined by 
changes in the polarisation angle, which T remains unaf- 
fected by. 

Since we are interested in using these results to con- 
struct a toy model of the full-sky synchrotron polarisation 
at high frequencies (30-100 GHz), there are two funda- 
mental questions that we need to ask. 

i) What is the effect of Faraday rotation on the Parkes 
data? 



As a result of propagation through an ionized medium, 
the direction of linear polarisation is rotated by an angle 
proportional to the radiation wavelength: <j> — RMA 2 . 

Gaensler et al. (2001) have made high resolution (1 
arcmin) polarimetric observation of a 28 deg 2 region in 
the Parkes field (325°. 5< I <332°.5, -0°.5< b <3°.5) at 
1.4 GHz. For this region, they derive a mean RM over the 
entire region of —12.9 ±0.1 rad m -2 . This at 2.4 GHz cor- 
responds to approximately 11°. However, the RM values 
span from 10 to 1000 rad m -2 , with typical uncertainties 
in RM measurements of ±20 — 40 rad m~ 2 . At a frequency 
of 2.4 GHz and resolution of 10.4 arcmin, a linear gradi- 
ent of 5 rad m _2 arcmin _1 in a foreground screen can sig- 
nificantly depolarise background radiation ( Sokoloff ct al 



1998J) and affect the measured variation of polarisation 



angle. Beam depolarisation is also significant at this fre- 
quency and resol ution for a le vel of dispersion in RM of 
30 rad nT 



le.g. 



Tribble 1991) 



Uaensler et al. (2UU1) conclude that the brightest po- 
larised features seen in this direction of the Galactic plane 
most likely represent intrinsic structure in the source of 
emission, but fainter structures are best explained as be- 
ing imposed by Faraday rotation on a uniformly polarised 
background by foreground material. 

This is likely to be the case also at 2.4 GHz and implies 
that a significant fraction of the angular variations which 
determine the power spectra of E and B for the Parkes 
data are likely to be induced by Faraday rotation. This 
scenario would predict that the angular power spectra of 
E, and B of synchrotron emission derived from the Parkes 
survey may be flatter than if they were derived from obser- 
vations at frequencies higher than 10 GHz, where Faraday 
effects are negligible (Beck 2001). 

ii) How representative are the results here derived on 
the angular power spectrum of linearly polarised emission 
at low Galactic latitudes of what may happen at high 
Galactic latitudes? 

Duncan et al. (1997) show that over much of the lon- 
gitude range covered by the survey, there exists a back- 
ground polarised component which appears to be inde- 
pendent of latitude (over the range of latitudes covered 



by the survey) . Besides the Parkes data, |Baccigalupi et al 



(2001)| also analysed five 10°-patchcs from th e Northern 
galactic plane survey by Duncan ct al. (1999)| at 2.7 GHz 



and three patches of different sizes at medium Galactic 
latitudes ( 2°. 5 < b < 17°. 5) from the survey of Uyamkcr 
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et al. (1999) at 1.4 GHz. By averaging the derived spec- 
tral indices for the polarised emission from all twenty fields 
they obtain a,L,z=ioo-800 = 1.8±0.3. They concluded that 
synchrotron polarisation maintains essentially the same 
statistical properties up to Galactic latitude b ~ 10°. 

However, by comparing high Galactic latitude fields 
with observations of the Galactic plane at 1.4 GHz Gray 



et al. ( -999) find that the magneto- ionic medium which in- 
duces the Faraday rotation on polarised emission appears 
to be concentrated in the disc component of the Galaxy 
(b ^ 10°). This would imply that higher latititude obser- 
vations of synchrotron polarised emission at low frequen- 
cies and for similar angular resolution may yield steeper 
angular power spectra. 

Indeed, if there is a variation of the slope with Galactic 
latitude one may expect this change to be in line with 
what is observed for the total emission, for which the 
angular power spectrum steepens with Galactic latitude 
Moreover, as pointed out by Davies fc Wilkinson (1999 



the magnetic field pattern may be more ordered at high 
Galactic latitudes, resulting in less small-scale structure. 

5. Full-sky synthetic maps of synchrotron 
polarisation 

As a starting point for creating a toy model of Galactic 
synchrotron polarisation at high frequencies we assume 
that the polarised component of synchrotron emission in 
a given direction is proportional to the unpolarised inten- 
sity, which in this case is expressed in terms of brightness 
temperature, T: 



Q = /Tcos(26>) 
U = fT sin(20) 



(1) 



where / is the fraction of polarised emission (or polar- 
isation degree) and 9 is the polarisation angle. As we 
have seen, at the frequency of Parkes observations this is 
clearly not a valid assumption, but it can be an acceptable 
starting point at frequencies higher than 10 GHz, where 
Faraday effects (of rotation and depolarisation) may be 
neglected. 

The brightness temperature spectrum of synchrotron 
emission is well described by a power law, 



T{v) oc v- 13 



(2) 



where (3 is referred to as the frequency spectral index. In 
the absence of non- uniform magnetic fields, the fraction 
of polarised radiation / is related to the spectral index 
( portiglioni fe Spoelstra 199% 

Therefore in order to generate synthetic full-sky polar- 
isation maps of synchrotron emission at a given frequency, 
we need at least three basic ingredients: a full-sky map of 
the total synchrotron emission at a "template" frequency, 
a full-sky map of the frequency spectral index and a model 



for the polarisation angle. The full-sky map of the spec- 
tral index will be used to generate / in a given direction 
and to extrapolate the total emission "template" to the 
other frequencies where Q and U are generated according 
to Eq. 1. 

5.1. Synchrotron total intensity (T) 



The 408-MHz survey of [Haslam et al. (1982)| provides, up 
to date, the best image of the total intensity of the full- 
sky Galactic synchrotron emission. At this low frequency 
the synchrotron emission dominates the Galactic diffuse 
signal ( Beuermann et al. 1985| ; Lawson et al. 1987). We 
used the 408-MHz survey processed by D.P. Finkbeiner, 
M. Davis and D. Schlegel (private communication). They 
have removed the point sources and de-striped the map by 
applying a Fourier filtering technique. Hereafter we will 
refer to this processed version of the 408-MHz survey as 
to the Cleaned 408 MHz map. 

The angular resolution of the 408-MHz survey is 0.85°. 
The angular resolution of Planck will vary from the 33 
arcmin FWHM in the 30 GHz channel to 5 arcmin FWHM 
in the 857 GHz channel. Since it is conservative to assume 
that synchrotron emission will present structure also at 
the smaller angular scale, we have superimposed artificial 
fluctuations at angular scale ^ 1° on the original map. 

The map of artificial fluctuations was obtained from 
a full-sky realisation of a Gaussian field having a power 
law spectrum with index a = 3 for I > 150 and vanishing 
to zero for I < 150. Before being added onto the Cleaned 
408 MHz^| the fluctuation map is weighted by the Cleaned 
408 MHz map itself} This operation can be viewed as an 
extrapolation of the spectrum of the Cleaned 408 MHz 
map beyond I — 150 with a power law spectrum with 
index a — 3. A power law spectrum with spectral index 
a = 3 is similar to that derived for the total synchrotron 
emission at high galactic latitude (|6| > 20° ) from the 
analysis of the three large sky radio survey up to I ~ 100 



(Bouchct fc Gispert 1999; Giardino et al. 2001) 



The Cleaned 408 MHz map with extrapolated fluctua- 
tions provides us with a realistic image of the full-sky total 
intensity of synchrotron emission at 408 MHz up to the 
degree angular scale and an estimate of how this emission 
may look like at the smaller unobserved angular scales. 
We use this map as the template for the Stokes T channel 
of the synthetic polarisation maps. The template at 408 
MHz is shown in Fig. |7[ 

3 convolved with a beam with FWHM of 0°.9 

4 A Gaussian realisation of a power law spectrum with index 
a — 3, for I > 150 and vanishing to zero for I < 150 has (by 
definition) power spectra with the same amplitude for all re- 
gions of the map. Such a map of artificial fluctuations (let us 
call it F) cannot be added onto the Cleaned 408 MHz (let us 
called this map H), where fluctuation amplitude for (/ < 150) 
is a strong function of the position of the sky. By multiplying 
F with H /ma,x(H) ("weighting") one obtains a map of fluctu- 
ations (at I > 150) for which the amplitude varies as a function 
of their position in a manner consistent with the Haslam map 
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Fig. 7. Mollweide projection of the Cleaned 408 MHz map with extrapolated fluctuations (i.e. the template for the 
Stokes parameter T of the synthetic polarisation maps). The map is in Galactic coordinates. The color scale has been 
histogram equalised. 



5.2. Synchrotron frequency spectral index 

Besides the 408-MHz survey the other two high quality, 
large sky radio surveys available to date are: the north- 
ern sky survey of [Reich fc Reich (1986)| at 1420 MHz 



and the more recent southern sky survey of Jonas et al 



(1998)| j The sky coverage of Lhe 1420-MHz survey and of 
the 2320-MHz survey are complementary Lu each oilier. 

The two surveys overlap over the Galactic declination 
range —20° < 6 < 32° in the right ascension inter- 
val 180° < RA < 360° and over Galactic declination 
range —20° < S < 13° in the right ascension interval 
0° < RA < 180°. Therefore we can combine the 408- 
MHz survey and the 1420-MHz survey in order to derive 
a map of the spectral index of radio emission in the north- 
ern sky and the 408 MHz and the 2326 MHz survey in 
order to derive a map of the spectral index in the south- 
ern hemisphere. We obtained three homogeneous maps in 
the following way: we took the surveys, as they are pub- 
licly available, resampled them into a HEALPix tessella- 
tion with nside= 256 and then median filtered with a box 
kernel of 9x9 pixels in order to suppress the point source 
signal. 

As we wish to determine the spectral index of the 
Galactic emission only, the extragalactic emission must 
first be subtracted. This is made up of two isotropic com- 
ponents, the CMB and the sum of unresolved extragalactic 
background. The CMB brightness temperature is 2.73 K. 
A phcnomcnological expression to estimate the integrated 
brightness temperature of extra galactic unresolved radio 



sources at different frequencies is given by Lawson et al 



408 MHz and in the 1420 MHz surveys. We have adopted 
it to evaluate this contribution at 2326 MHz. The total 
amount of extragalactic contribution for the three surveys 
together with basic information about their uncertainties 
is summarised in Table |[ For the 1420 M Hz survey we 
adop ted a baseline correction of —0.13 K ( Lawson et al 



1987) 



After subtraction of the extragalactic contribution, the 
brightness temperature frequency spectral index f3 v i/ U 2 of 
the Galactic radiation can be calculated. It is simply given 
by: 



T Ul /T U2 
log(^i/i/ 2 ) 



(4) 



1987) . Lawson uses it to estimate this contribution in the 



Over the whole observed region the average value of 
the spectral index obtained by combining the 408 MHz 
and the 1420-MHz survey is (/34os/i42o) — 

2.78 ± 0.17. 

By combing the 408 MHz and the 2326-MHz survey 
one obtains an average value of the spectral index of 
(/3408/2326) = 2.75 ± 0.12, over the whole region observed. 
Here, the error is only the formal rms of the spectral index 
maps and does not account for the systematic uncertain- 
ties. The mean values for (3 agree well with the value found 
by [Platania et al. (1998)| of (/3 40 8/380o} = 2.73 ± 0.08. 

All the surveys are affected by scan-to-scan baseline er- 
rors. When the ratio of brightness temperatures at two dif- 
ferent frequencies is calculated the effect of these residual 
errors is enhanced. Because of this, the map of /?4os/i420 
and the map of A108/2326 d° not match in the overlapping 
region. We tried various techniques to look for a base- 
line offset for one of the three input maps that would 
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Survey 


Scale 


Base level 


Extragalactic 


Baseline 


Reference 


Frequency (MHz) 


error (per cent) 


error (K) 


background (K) 


correction (K) 




408 


5 


3 


5.92 




Lawson et al. (1987) 


1420 


10 


0.6 


2.83 


-0.13 


Lawson et al. (1987) 


2326 


5 


0.080 


2.75 




Jonas et al. (1998) 



Table 2. Characteristics of the three radio continuum surveys that we used to produce a full-sky map of the radio 
brightness temperature spectral index. 



minimise the differences between the two spectral index 
maps. However discontinuities between the two spectral 
index maps always remained apparent. 

In order to construct the full-sky map of the spectral 
index, therefore, we used A108/1420 f° r regions of the sky at 
declination 6 > 13° and /3408 / 23 26 f° r regions at declina- 
tion 5 < —20°. In the overlapping region of the 1420 and 
2326-MHz surveys (-20° < S < 13°) we set the spectral 
index to: 

/^overlap = 5/?408/1420 + ^/3408/2326 

where g = exp[-0.5(<5 - 13°) 2 /(10°) 2 ] and h = 1 - g. 
This choice of the weights g and h ensures that the two 
/3-maps merge smoothly into each other. The region of sky 
at declination S < 83° remains un-observed at 2326 MHz. 
So we "padded" this region with the spectral index value 

(Al08/2326)- 

The spectral index map obtained in this way is how- 
ever still affected by the striation due to the scan-to-scan 
baseline errors in the input data. It is possible to remove 
the effect by Fourier filtering the data but, because the 
scanning direction differs from one survey to the other, 
one would not then be able to combine like with like. 



Following Lawson et al. (1987) we adopted the simpler 
procedure of convolving the data to a lower resolution, so 
that the baseline errors are averaged sufficiently for the 
striation on the color-scale to disappear. The final map of 
the spectral index is at a resolution of 10° (at FWHM), 
as shown in Fig. ^. 

5.3. Synthetic polarisation angle 

We found that the power spectra of the E-B transform 
of the sinus and cosinus of the polarisation angle of the 
Parkes data are described by power laws with spectral in- 
dices of ct e = 1.74 ± 0.14 and a b = 1.69 ± 0.12 in the 
multipole range I = 40-250. One can construct a full-sky 
map of a random angle having E-B transform of its sinus 
and cosinus with the same spectra of the e and b compo- 
nent of the Parkes data in the following way. Generate two 
full-sky realisations of a Gaussian field characterised by a 
power law spectrum with spectral index a = 1.7. Let us 
call these maps x and y. Normalise x and y to the interval 
[—1,1]. Generate the full-sky map of the random angle by 
computing 9 = \ arctan(x/n, y/n), where n = \J x 2 + y 2 . 
We used the random angle map obtained in this way as 
the polarisation angle of our toy model of Galactic syn- 
chrotron emission. 



Fig. |9| shows the map of the polarisation angle derived 
from the Parkes data and the map of the synthetic polar- 
isation angle for the same region of sky. The structure of 
the angle fluctuations in the two maps looks similar. 

5.4. The synthetic maps at 30 and 100 GHz 

In order to generate the synthetic maps of synchrotron 
emission and linear polarisation at higher frequencies the 
total intensity map at 408 MHz has to be extrapolated. 
One can use the spectral index map directly, assuming a 
fixed power law to higher frequencies, or use it in con- 
junction with the local electron spectrum, as in Bennett 



t al. (1992) , in order to derive maps of (3 between any two 



given frequencies and to attempt to allow for the observed 



steepening of (3 with increasing fr equency ( Lawson et al 
19871; [Banday fc Wolfendale 1990|). 



For the local electron spectrum, we used the poly- 
nomial fit to the spectral shape of electron measure- 

for energies 0.1 < 



Bennett et al. (1992) 



ments given m 
E(GeV) < 100, and likewise a power law E~ 3 312 dE 
for E > 100 GeV. By interpreting the spatial variation 
in (3 (408-1420 MHz and 408-2326 MHz) as due solely 
to variations in a synchrotron effective Galactic mag- 
netic field, -B e ff) we derived a full-sky map^| of B c g , with 
(_B c ff) = 2.0 ± 2.3 /iG and a median value of 1.4 /iG. This 



is remarkably close to the values derived by Bennett et al 
|(1992)| who did not have the data from the 2326-MHz sur- 
vey for the southern hemisphere. 

The same expression for the local electron spectrum 
and the map of synchrotron effective Galactic magnetic 
field were then used to obtain the maps of the spectral 
index between 408 MHz and 30 and 100 GHz, necessary for 
extrapolating the total emission map. The (3 maps mean 
values are (/3 .408/3o) = 2.91 ± 0.09 and (/3 .408/ioo> = 
2.96 ± 0.08. The value of (A).408/3o) is consistent with the 
upper limit of (3 > 3 between 408 MHz and 31.5 GHz 
obtained by Kogut et al. (1996) and the valu e of (3 ~ 2.8 
between 1420 MHz and 19 GHz derived by |de Oliveira 
|Costa et al. (1998)| . 

The Q and U channels of the synthetic linear polar- 
isation maps of Galactic synchrotron emission at 30 and 
100 GHz were then obtained by combining the intensity 
map (T) extrapolated at 30 and 100 GHz with the corre- 
sponding spectral index map and the random angle map, 



5 nearly full-sky: only the region at declination S < 83° re- 
mains un-observed at 2326 MHz 
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Fig. 8. Mollweide projection of the full-sky map of radio spectral index obtained by combining the 408-MHz survey, 
the 1420-MHz survey and the 2326-MHz survey (see text). The map is in Galactic coordinates. 
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Fig. 9. Polarisation angle from the Parkes survey (top) and from the synthetic map of synchrotron polarisation 
(bottom). Maps are in Galactic coordinates. 



The Q and U channels 
The three input maps 



according to Eq. (1) and Eq. 
at 30 GHz are shown in Fig. 
are sampled into a HEALPix tessellation with a pixel size 
of 6.9 arcmin (nside = 512). Therefore also the synthetic 
maps of T, Q and U are at a pixel resolution of 6.9 arcmin. 



Because our toy model for / is a very slowly varying 
function of the position in the sky the maps of Q and U 
derived according to Eq. (1) have temperature fluctuations 
mostly determined by the polarisation angle. Therefore, 
like for the Parkes data, the angular power spectra of the 
E and B components of the synthetic linear polarisation 
maps reflect the spectra of e and b. They are power-law 
spectra with indices «_b = 1.67±0.09and as = 1.73±0.09 
over the I range 40-250 (ag = 1.59 ± 0.03 and olb = 
1.68 ± 0.03 for the I range 3-250, that can be probed since 
in this case the power spectra are derived from full sky 
maps). These values are very close to the values derived 



for the spectra of the E and B components of the Parkes 
data. 



5.4.1. Comparison of the power spectra with 
predictions for the CMB 

For CMB studies one is particularly interested in the prop- 
erties of radio emission at high Galactic latitudes because 
this is where the Galactic signal is weaker and where sen- 
sitive CMB measurements can be made. Therefore, we 
computed the angular power spectra of the T, E and 
B components of the extrapolated maps with a Galactic 
plane cut-off at latitude \b\ = 20°. The results are shown 
in Fig. [H], where they can be compared with the sig- 
nal expected from the CMB. The CMB spectrum has 
been computed assuming a fiat inflationary model with 
purely scalar scale- invariant fluctuations, vacuum energy 
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Fig. 10. The Q and U channel of the synthetic map of Galactic synchrotron linear polarisation at 30 GHz. The map 
is in Galactic coordinates. 



Oa = 0.70, cold dark matter density /i 2 f2 c d m 



0.12 and 



baryon density h = 0.024, (where h = 0.7 is the as 



sumed Hubble parameter) , as derived by Netterfield et al 
(2001)|from combining the latest CMB observations with 
constraints derived from measurements of the large scale 
structure and results from recent measurements of type la 
supernovae. 

The top panel of Fig. [ll] shows the angular power spec- 
tra of the T components of the two extrapolated maps and 
of the CMB. Note the absence of a perceivable change in 
the slope of the power spectra of the T components of the 



synthetic maps at I = 150, the multipole order at which 
artificial signal was added to the observed signal at 408 
MHzQ. 

From the figure, it is apparent that at 100 GHz the con- 
tribution to the total sky signal from high latitude syn- 
chrotron emission is negligible. At 30 GHz synchrotron 



6 In order to prevent the aliasing noise that effects power 
spectra with spectral index J; 3 when a Galactic cut is applied 
we used an apodysed Galactic cut off and the modified power 
spectra definition given in Giardino et al. 2001, in which the 
terms <2;o are set to zero 
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emission can be a significant contribution to the sky 
anisotropy at scales larger than ~ 20°, but it is more 
than an order of magnitude weaker than the computed 
cosmological signal at degree angular scales where the 
first "acoustic" peak of the CMB is observed. This con- 
firms the results discussed in 



Giardino et al. (2001) 



that 

were obtained with a different approach using the index 
of the angular power spectra of the survey at 2326 MHz 
by ponas et al. (1998)| and the COBE DMR upper limit 
to synchrotron temperature fluctuations at 31.5 GHz at 



the 7° scale. The COBE DMR upper limit ( [Kogut et al 
1996(35 indicated in the figure by a pentagonal point. 

The angular power spectra of the E and B compo- 
nents of the synthetic maps are shown in the lower panel 
of Fig. O. At small angular scale the power of temperature 
fluctuations in the E and B components are mostly de- 
termined by change in the polarisation angle rather than 
a change in the degree of polarisation. This is the reason 
why temperature fluctuations in the E and B components 
can become more intense than temperature fluctuations in 
the T channel. 

The figure shows that at 30 GHz the high latitude syn- 
chrotron signal may dominate the cosmological polarised 
signal over the entire multipole range. The extrapolated 
signal at 30 GHz is consistent with the current upper limit 
on linearly polarised signal of 10 /iK at the 7° scale in the 
26-36 GHz band (Keating et al. 2001) and with the up- 
per limit of 16 /iK at degree scale in the 26-46 GHz band 
( [Nctterficld et al. 1995| ). These upper limits also refer to 
regions of the sky at galactic latitudes \b\ > 20°. 

At 100 GHz the synchrotron signal can make a signifi- 
cant contribution to the sky polarisation at large angular 
scale but it is significantly weaker than the CMB E com- 
ponent at sub-degree angular scale. Therefore even using 
the most conservative assumptions one can conclude that 
the contribution of synchrotron polarisation to the sky sig- 
nal will not hinder detection of the CMB E component, if 
this is present at sub-degree angular scale at the intensity 
predicted by current cosmological models. 

From the diagram the absence of the beam cut off in 
the spectra of the synthetic maps is noticeable. These are 
artificial maps and therefore they were generated at full 
resolution. 



5.5. Scope and limitations of the synthetic maps 

The synthetic full-sky maps of the Stokes parameters T, 
Q and U of Galactic synchrotron emission at 30 and 100 
GHz are a first attempt to model the full-sky synchrotron 
polarisation at high frequencies. The maps are needed as 
input for the simulation activity that is carried out in 
preparation for the Planck missions. 

To study the effects of polarised foreground on MAP 
observations of the CMB, synthetic maps of polarised 
synchrotron em ission have been generated by Kogut fc 
Hinshaw (2000) . They also assumed the polarised signal 
to be proportional to the total intensity T, with the Stokes 
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Fig. 11. The power spectra of the T, E and B compo- 
nents of the synthetic maps of Galactic synchrotron po- 
larisation at 30 and 100 GHz (dashed lines), for |6| > 20°, 
compared with the signal expected from the CMB (con- 
tinuous line). The point with pentagonal shape indicates 
the COBE DMR upper limit to synchrotron temperature 
fluctuations at 31.5 GHz. The horizontal line is the up- 
per limit to polarised signal by Keating et.al (2001) in the 
frequency band 26—36 GHz. The triangular point is the 
upper limit to polarised signal by Netterfield et al. (1995) 
in the frequency band 26—46 GHz. 



Q and U components given by Eq. 2. For the polarisation 
degree, /, they assumed a varying Gaussian field such that 
(/) =0.1. They constructed the polarisation angle start- 
ing from Gaussian random fields, x and y, with flat power 
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spectra^], which lead to E and B components whose spec- 
tra have indices of a ~ 0. The polarisation maps were only 
generated for a fiducial frequency of 40 GHz. 

We have constructed the synthetic maps by combin- 
ing in a conservative way many of the observational data 
currently available on synchrotron total emission and po- 
larisation. The maps that can be obtained at higher fre- 
quencies by extrapolating the total intensity map with the 
spectral index map are consistent with available observa- 
tional upper limits on synchrotron emission and polari- 
sation at frequencies higher than 10 GHz. There remain 
however some inconsistencies between the synthetic po- 
larisation maps extrapolated at 2.4 GHz and the Parkes 
data. 

The synthetic polarisation maps reproduce the slope 
of the angular power spectra of the E and B components 
of the Parkes survey, but not their amplitude at 2.4 GHz. 
Moreover, the power spectra of T, E xT and L differ both 
in shape and normalisation to the ones of the Parkes data. 
This is because of the way that the toy model has been 
constructed. We derived information on the polarisation 
angle (and the E and B power spectra) from the Parkes 
survey, while for the intensity channel T we used the data 
from the Haslam survey and for the polarisation fraction / 
we assumed the theoretical value given by the synchrotron 
spectral index in a uniform magnetic field. The reasons 
behind these choices are the following. 

As already mentioned, the Haslam map provides the 
best image of the total intensity of the full-sky synchrotron 
emission, to date. It is natural to adopt it as a template 
for T as opposed to using the T channel of the Parkes 
survey. This is limited in sky coverage and being at low 
Galactic latitude and a higher frequency contains a higher 
fraction of thermal emission from ionised gas. 

The Parkes survey provides us with information about 
the statistical properties of the polarisation angle, so we 
adopt it for this purpose. Using the statistical properties 
of the polarisation angle derived from a survey at deci- 
metric frequency is a conservative approach (i.e. it will 
unlikely lead to an underestimate of the level of fluctu- 
ation of polarised synchrotron emission). This is because 
Faraday rotation of a synchrotron background due to a 
foreground "screen" of ionised medium will tend to in- 
crease the level of fluctuations in the observed polarisation 
direction. On the other hand Faraday differential rotation 
may decrease the polarised fraction of the synchrotron sig- 
nal at the Parkes frequency. Therefore using the amplitude 
of the E and B components of the Parkes data may lead 
to an underestimate of the polarisation fraction at higher 
frequencies (Jb 10 GHz). 

Calculating the theoretical fraction of synchrotron po- 
larisation by assuming a completely uniform magnetic 
fields could be considered an overly conservative approach 



because it is known that the interstellar medium is per- 
vaded by an isotropic random ("turbulent") field. In the 
presence of a random magnetic field the observed level of 
polarisation is given by / D b s — d - f. Here / is the intrinsic 
level of polarisation (given by Eq. [|) and d is a coeffi- 
cient (< 1) which depends on the relative intensity of the 
random field, -Bturb, and the field which is regular in the 
whole resolution element, i? rcg (Burn 1966; Beck 1998). 



7 the fields were then convolved with a Gaussian with vary- 
ing FWHM, as the authors were interested in investigating the 
effect on CMB observations of the eventual presence of a co- 
herence angle in polarised Galactic emission 



Synchrotron polarisation observations at decimeter 
wavelength imply a ratio of regular to total field strengths 
within abou t a kpc from the Sun of (-B reg / -Bturb} — 0.6 
( |Beck 200l|) , which, in the popular assumption of energy 
equipartition between cosmic rays and magnetic fields 
would imply an average level of depolarisation of (d) = 0.3. 
However the actual value varies substantially between 
points and may depend on the size of the resolution ele- 
ment. These values are derived from decimeter wavelength 
surveys with degree-scale resolution. The level of depolar- 
isation by the random component of the magnetic field 
can be expected to decrease as the sky is observed at 
higher resolution and the projected size of the beam ap- 
proaches the field coherence length. Indeed convolution of 
our synthetic maps with a Gaussian beam with a FWHM 
of 1° yields an average level of depolarisation of (d) = 0.5. 



6. Summary and Conclusions 

In order to obtain information about the statistical prop- 
erties of synchrotron polarised emission, we have derived 
the global angular power spectra of the Parkes survey, a 
radio continuum and polarisation survey of the Southern 
galactic plane at 2.4 GHz. The angular power spectrum 
of the polarised intensity (L) is well approximated by a 
power law with fitted spectral index = 2.37 ± 0.21 in 
the multipole range I = 40-250 that can be probed given 
the Parkes survey sky coverage and angular resolution. We 
exclude that up to multipole order I = 150, the spectral 
index value is affected by the presence of discrete signals 
such us point sources or instrumental noise at pixel scale. 

The angular power spectra of the E and B components 
of the polarised emission is significantly flatter than the 
spectrum of L. They have fitted spectral indices respec- 
tively of a E = 1.57 ± 0.12 and a B = 1.45 ± 0.12, in the 
multipole range I — 40-250. This is because temperature 
variations in the Stokes channels Q and U of the data 
are mostly determined by variation in polarisation angle 
rather than polarisation intensity. The E-B transform of 
the sinus and cosinus of the polarisation angle, the e and 
b components, are in fact also well described by power 
laws with fitted spectral indices of a e = 1.74 ± 0.14 and 
on, = 1.69 ±0.13, in the same multipole range I = 40-250. 

We generated a full-sky map of a random polarisation 
angle having e and b components with power spectra with 
the same spectral indices of the e and b components of 
the Parkes data. We combined this map with a template 
of synchrotron emission at 408 MHz, derived from the 408- 
MHz survey, and with a full-sky spectral index map, that 
we obtained by combining the 408-MHz survey with the 
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1420-MHz survey and the 2326-MHz survey. The simple 
recipe on how to combine the three maps in order to derive 
a map of linear polarisation is given by Eq. (1) and Eq. (||). 

In this way we have constructed synthetic maps of the 
Stokes parameters (T, Q and U) of Galactic synchrotron 
emission at 30 and 100 GHz. The angular power spec- 
tra of the E and B components of these synthetic maps 
have slope very similar to the E and B components of the 
Parkes data. Their polarisation intensity is proportional 
to the total intensity T as determined by the frequency 
spectral index in the conservative hypothesis of a uniform 
magnetic field. 

We have compared the angular power spactra of the 
synthetic maps with predictions for the CMB and con- 
clude that while direct observations of the CMB E com- 
ponent at 30 GHz may be precluded by synchrotron polar- 
isation, at 100 GHz, this component of Galactic emission 
will not hinder the observation of cosmological polarisa- 
tion, if this is present at the predicted levels. 

The synthetic maps of the linear polarisation of 
Galactic synchrotron can be useful as a toy model to study 
the effect of the polarised foreground on planned observa- 
tions of the CMB. Moreover they can serve as the basis for 
the development of more advanced models which incorpo- 
rate data on the structure of the Galactic magnetic field. 
In turn these refined models will be of fundamental im- 
portance in interpreting the microwave polarisation data 
that will be provided by the Planck and MAP satellite 
missions and the SPORT experiment. 

The synthetic maps of the Stokes parameters 
T, Q and U of Galactic synchrotron emission at 
30 and 100 GHz together with the full-sky spec- 
tral index map are available for down loading at 
ftp : / / astro . esa . int/pub/ synchrotron. 
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